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Abstract 

We analyze CP violation in the supersymmetric standard model (MSSM) 
embedded minimally into a left-right symmetric gauge structure with the 
seesaw mechanism for neutrino masses. With the plausible assumption of 
universal scalar masses it is shown that CP violation in the hadron sector of 
the MSSM is described by a single phase residing in the supersymmetry break- 
ing Lagrangian. This improves the CP properties of the MSSM by providing 
a natural solution to the SUSY CP problem. Furthermore, e'/e vanishes in 
this model above the seesaw scale; extrapolation to the weak scale then leads 
to a prediction in agreement with the NA31 and the recent KTeV observa- 
tions. The electric dipole moment of the neutron is naturally suppressed, we 
estimate its magnitude to be about 4 x 10 -29 ecm. Additional predictions 
include a tightly constrained super particle spectrum and vanishingly small 
CP asymmetries in the B meson system. 



Recent evidences for neutrino oscillations imply that the standard model must be ex- 
tended to accommodate small neutrino masses. An elegant model that provides an expla- 
nation of the small neutrino masses via the seesaw mechanism |J is the left-right symmetric 
model of weak interaction. Two seemingly unrelated puzzles of the standard model, viz., 
the stability of the Higgs-boson mass and the question of the origin of the electroweak sym- 
metry breaking, seem to require its supersymmetric extension - the MSSM - for a proper 
resolution. It is well known that MSSM by itself is plagued with a number of problems 
which were not present in the standard model. In this Letter, we will be concerned with one 
such problem, viz., the emergence of a plethora of arbitrary CP phases that require severe 
fine-tuning to explain the smallness of observed CP violation in the Kaon system as well as 
the non-observation of the electric dipole moment of the neutron |]3|-|| • This is known as the 
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SUSY CP problem. A promising approach towards its resolution would be to realize MSSM 
as a low energy limit of a theory where this as well as other problems of MSSM are cured. 
Here we explore supersymmetric left-right model with the seesaw mechanism for neutrino 
masses as a candidate theory. 

It has been pointed out that in the supersymmetric version of the left-right model (SU- 
SYLR) |||7j which embodies the seesaw mechanism (and therefore explains the small neu- 
trino masses) the constraints of parity invariance provide a simple resolution of the SUSY CP 
problem ||. Neutrino masses as suggested by SuperKamiokande require that the scale above 
which the SUSYLR model manifests itself is at least 10 12 GeV, if the seesaw mechanism is 
implemented via the renormalizable terms in the superpotential, or close to the canonical 
GUT scale of 10 16 GeV if one envisions the Planck scale suppressed non-renormalizable 
terms as the source of the seesaw. In either case the model could finally be embedded into 
an SO(10) grand unified theory. The basic assumption of this paper is that above the scale 
where the SUSYLR model is valid, the Lagrangian is invariant under the parity transforma- 
tion and below the seesaw scale the theory is MSSM, with constraints on its parameters as 
required for its embedding into the SUSYLR model. 

We will focus our attention on the flavor mixing and CP violation in the minimal versions 
of the high scale SUSYLR (or SO(10)) model. By minimal we mean that we must have only 
one multiplet that gives rise to fermion masses, i.e. one left-right bidoublet (10 in the case 
of SO(10)). Since this one multiplet contains both the H u and multiplets of the MSSM, 
it immediately leads to a proportionality between the up and the down Yukawa coupling 
matrices P,|10|. This is called up-down unification. The quark mixing angles all vanish at the 
tree-level, but they are induced by loop diagrams involving the exchange of supersymmetric 
particles. This considerably restricts the flavor and CP violating interactions in the model 
and makes it very predictive. Note that we have not assumed the existence of any extraneous 
discrete symmetries. This is one of our key starting points. 

There is a second set of constraints on the model which follows from the existence of 
parity invariance of the Lagrangian prior to symmetry breaking || . Above the seesaw scale 
they imply that the familiar fi and B parameters as well as the gluino mass terms are 
real. Furthermore, the Yukawa coupling matrices as well as the SUSY breaking A matrices 
(the trilinear terms that involves the squarks) are Hermitean. In fact, left-right symmetry 
implies that there is only one A matrix in the squark sector, which evolves to the two A u ^ 
matrices of low energy MSSM. It was noted in Ref. |§ that these parity-implied constraints 
solve the SUSY CP problem in the sense that the electric dipole moment of the neutron 
is comfortably consistent with the present experimental limits Jll| without any need for 
fine-tuning. We thus see that the combination of up-down unification and the constraints 
of parity invariance considerably restricts flavor structure of both the quark and the squark 
sector of the model PJlOl and one might therefore expect that in addition to solving the 
SUSY CP problem, there are predictions by which the model can be tested. 

In this paper we supplement the seesaw-constrained MSSM just described with the 
plausible assumption of universal scalar masses. We shall keep the trilinear scalar A terms 
arbitrary, subject of course to left-right symmetric constraints. In this theory, there is only 
one CP phase residing in the A term that characterizes all CP violating phenomena in the 
quark sector. It is thus on par with the conventional CKM model as far as the number of 
CP violating phases is concerned. 
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At tree level, the up and down Yukawa matrices in the model are proportional. This 
results in the vanishing of the CKM angles at tree level. Non-zero quark mixings arise only 
from one-loop corrections involving the elements of the A matrices. There is no other source 
of flavor mixing in the model. The are thereby fixed to a narrow range, resulting in a 
very predictive model. The value of the single CP phase is fixed by the requirement that the 
model reproduce the observed value of ex, the indirect CP violating measure in the Kaon 
system. The resulting low energy theory is the MSSM, but without the SUSY CP problem 
and with its parameters restricted to a very narrow range. 

An interesting prediction of the model is that, owing to the constraints of parity invari- 
ance, e'/e is vanishingly small above the seesaw scale vr. However, as the theory is RGE 
evolved to the weak scale, manifest parity invariance disappears and a non-negligible value 
for e' / e emerges. We calculate this value and find it to be in good agreement with the recent 



KTeV [|T|] and previous NA31 |13| results. 

A similar suppression also occurs for the electric dipole moment (edm) of the neutron and 
we find that its value at the weak scale is ~ 4 x 10~~ 29 ecm. This and the fact that up-down 
unification restricts the allowed parameter space of the MSSM considerably provide tests of 
the model. 

Let us start by giving a brief derivation of the up-down unification relation in the SU- 
SYLR models. As is well known, the gauge group for this model is SU{2) L x SU{2) R x 
U{1)b-l with the standard assignment where Q,Q C denote left-handed and right-handed 
quark doublets and $ denotes the (2,2,0) Higgs bi-doublet. The SU(2)r x U{1)b-l sym- 
metry could either be broken by B — L = 2 triplets - the left-handed triplet A and the 
right-handed triplet A c (accompanied by A and A c fields, their conjugates to cancel anoma- 
lies) or by B — L = 1 doublets x (left) and x c (right) along with x an d X c - Let us write 
down the gauge invariant matter part of the superpotential involving these fields: 

W = Y g Q T T 2 $T 2 Q c + Y 1 L t t 2 ^t 2 L c 

+ i(iL T r 2 AL + f c L cT r 2 A c L c ) . (1) 

Below the vr scale, the B u ,d contained in the bi-doublet field will emerge as the MSSM 
doublets, but in general in arbitrary combinations with other doublet fields in the model. 
The single coupling matrix Y q therefore describes the flavor mixing in the MSSM in both 
the up and the down sectors leading to the relations 

Y u = lY d ; Y t = 7 Y v u (2) 

which we call up-down unification. The parameter 7 is unity if the multiplets H u and Hd of 
MSSM are contained entirely in $(2,2,0), but 7 can differ from one if additional doublets 
contribute to H u and Hd- At first sight the first of the relations in Eq. (2) might appear 
phenomenologically disastrous since it leads to vanishing quark mixings and unacceptable 
quark mass relations. We showed in Ref. |9j that after including the one-loop corrections 
involving the exchange of supersymmetric particles to the relations in Eq. (1), there exists 
a large range of parameters (though not the entire range in the MSSM) where correct quark 
mixings as well as masses can be obtained. In Ref. |J, we explored the parameter space 
that allowed for arbitrary bilinear squark masses and mixings as well as arbitrary form for 
the supersymmetry breaking trilinear A matrix. We focused on a class of solutions for large 
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tan/3 ~ 35 — 40, corresponding to 7 =1, where all quark masses mixings and CP violating 
phenomena could be explained. The magnitude of tan/3 can be reduced for 7 > 1. In this 
paper, we focus on a predictive scenario where all flavor mixing arises from the trilinear A 
terms. We use the small tan/3 scenario (7 ^> 1) to explain the observed CP violation while 
satisfying the FCNC constraints from the K meson system. 

The model becomes much more predictive once the assumption of universal scalar masses 
is imposed on the theory. This implies that at the Planck scale, the only phase of the SU- 
SYLR theory resides in one of the three off-diagonal entries of the SUSY breaking trilinear 
coupling matrix A. To see how this comes about, note that parity symmetry of the La- 
grangian imposed above the vr scale implies that Arg(B) = Arg(fi) = and Y Q) i = 
and Agj — A' ql . A Hermitean 3x3 matrix has only three independent phases. However, 
above the vr scale, quark masses can be diagonalized and made real (this is because there 
is only one matrix Y q in Eq. (1).) It therefore follows that redefining the phases of any two 
quark superfields (in both the right and the left sector) we can make two of the three off 
diagonal A^-'s real. Thus we are left with only one phase in the theory. Thus the tree level 
MSSM parameters above the vr scale can be summarized as follows (we only discuss the 
quark-squark sector): 

AfW = 7 tan/3Mi 0) (3) 
M| = Ml = M} c = Ml ud = mil; 
M w = Mg = M~ g = Mi/2 

f A n A 12 A l3 e iS ^ 

A 12 A 22 A 23 I . (4) 
\A 13 e i5 ^ A 23 A 33 



A 



The parameter 7 is related to the mixing between the SU(2)i doublets in $ and those in 
other multiplets in the theory such as \ in the case of the non-renormalizable seesaw model 
or those in (2, 2, ±2) multiplets which may be included in the renormalizable seesaw modelQ. 
Note that all A%j in Eq. (4) are real and further, we could have chosen to place the phase 

at any off diagonal entry of A. However the final results are independent of this choice. 

The first task before us is to compute the one-loop corrections to the quark mass matrices 
both in the up and the down sector and obtain the desired masses and mixings. The one-loop 
expressions for the mass corrections to the down type quarks are given in the Refs. P, |i~0| , |T4"|| . 
The up sector corrections can be obtained by replacing (Ad, Xd, Vd) by (A u , X u , v u ). In 
the down sector, there are three types of flavor contributions: Md = Vd[Y q (l + citan/3) + 
C2(Ad/ Msusy) + £3^33]. Here the q are dimensionless loop factors. The c\ and c 2 terms 
arise from the gluino graph, the c 3 term which contributes significantly only to the 6-quark 
mass is from the chargino graph. M u is given bjM u = v u [Y q (l + cij tan (3 + c 2 (A u / Msusy)]- 
Clearly there is a mismatch between M u and M^, which implies violation of proportionality 
and non-zero CKM angles. 



x The advantage of including the (2, 2, ±2) multiplets instead of the x _ type doublets is that the 
former maintain the property of automatic R-parity conservation. 
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Although the scalar masses are assumed to be universal at the Planck scale, the non- 
diagonal nature of the A matrix will induce via the RGE off-diagonal elements in the up and 
down squark mass matrices. Since our calculation is going to be done at the SUSY scale of a 
few hundred GeV, we must extrapolate the parameters down from the Planck scale via the 
vr scale down to Msusy- The RGE's for extrapolation below Vr are those of MSSM and are 
well known ||15|| . Between vr < fi < Mpg, we use the RGE corresponding to the SUSYLR 



model. We keep only the one-loop terms. It turns out that the main effect of running from 
the Planck scale to v r in the squark sector is to split the third generation squarks slightly 
from the first two generations due to the large third generation Yukawa coupling. This 
effect is further amplified via the RGE in the process of running from vr to the SUSY scale. 
We work in a basis where the tree-level Yukawa coupling matrix is diagonal. As a result, 
Yukawa matrices at Msusy also will be diagonal. However, the superpartner masses which 
start at the Planck scale as diagonal matrices acquire off-diagonal terms both at the vr scale 
and at M S usy- At M SUSY , the Q, u c and d c masses are no longer equal. The off-diagonal 
entries of these matrices become complex. Similarly, the A-matrix which was Hermitean at 
the vr scale also loses its hermiticity. Using these extrapolated quantities, we compute the 
one-loop corrections to the up and down quark mass matrices and diagonalize them to fit 
the known quark mixings and masses. We find a range of input parameters at Mpi which 
leads to the correct quark masses and mixings. We then make sure that the chosen values 
of the Aij's do not lead to excessive flavor changing neutral current effects. We have used 
the constraints quoted in Ref. |16| , |17| . 

Before presenting an explicit example of the parameter choice and the corresponding 
mass predictions, let us recall one of the main results of Ref. 0. We showed that, to get 
the correct mixing angles, we need to have 5i2,lr — 4.4 x 1CT 3 — 6.2 x 1CT 3 , 5 2 3,lr — 
(.84 — 1.8) x 1CU 2 , where the range comes from varying the parameter x = M?/m?. 8\ 2 ,lr is 
determined from the Cabibbo angle, 523,lr from V c b. (Here Sij are the flavor violating squark 
mixing parameters.) On the other hand, from Ref. [jnj, we note that the upper bound on 
Itji{8\2,lr) — (2 — 4) x 10 -4 . It is then clear from this that the phase 5f 3 in the A matrix 
should be of order 10~ 2 . Our detailed numerical analysis also seems to generate fits to all 
parameters only for a phase of this order of magnitude. Furthermore, we find that if we scale 
the squark masses, the A matrix and Mg by a common factor k, the quark mixings remain 
unchanged. This might lead one to suspect that the SUSY breaking parameter range is not 
limited in the theory. But since ex in our model is coming entirely from the SUSY box graph 
and it scales like m^ 2 , the squarks cannot be too heavy. There is a scaling relation between 
the CP phase and the SUSY breaking parameters. The e'/ e however scales differently. As a 
result, we are forced to a narrow range of the SUSY breaking masses. 

To present a concrete example, we consider a case with tan/3 = 3, m = 80 GeV, 
and Mi/2 = 180 GeV. For the trilinear A matrix we choose: (An, Als, A22A23, A 33 ) = 
(1.2, 1.8, —2.2, —12, 17, 50) GeV and 5f 3 = 0.02. We determine /1 from the radiative breaking 
of the electroweak symmetry. Its magnitude in the above parameter space is fi ~ 290 GeV. 
We choose the sign of fi as preferred by b — > sj decay. For the quark masses and mixings, 
we find: 

V us ~ -0.21, V cb ~ 0.035, V ub ~ -0.0033, V td ~ -0.012, and J-~7x 10~ 7 . (5) 
M d = (-0.0042,-0.059,2.62) GeV, M u = (-0.0024,0.61,162) GeV. 
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The top mass (pole) is 172 GeV. The other masses at their respective mass scales (or at 1 
GeV for u, d, s) can be obtained by mulityplying with the following QCD correction factors 
rj: 77&=1.59, r/ c =2.1, i] s =2A, 77^=2.4 and r] u =2A. The fit for the quark masses and mixings 
is quite satisfactory. 

The squark mass matrices are 6x6 with the 3x3 submatrices denoted by Mf jL ,M RR and 
M\ R . These submatrices need to be rotated by the matrices which are used to diagonalize the 
up and the down quark mass matrices. If the left-handed and the right-handed rotations for 
the quark masses are given by Ui » and U r i where i's can be u or d, then we write the rotated 



down type squark submatrices (for example) as: Ui^M^Uj d , U r ^M RR Ul d and \J\^M\ R \j] 



r,d' 



We present only the rotated down squark mass matrices in this case since they are the only 
ones that enter in the discussion of CP violation in K-decays: 



M 



LL 



M 2 RR 



M LR 



I 210442 -1.8 - 198* 

-1.8 + 198i 209509 
V -3.5 + 398?: -1594 + 0.86?: 

/ 193030 -2.4 - 263.6?: 

-2.4 + 263.6?: 191835 
V -5.7 + 651?; -2583 + 1.4? 

/ 180 + 7.2 x 10" 
0.23 -86.3?; 
7.9 - 874.8? 



-3.5 - 398? 
-1594- 0.86?; 
174085 

-5.7-651? 
-2583 - 1.4? 
187244 

0.23 + 86.2?; 8.85 + 979.8?; 
2430 + 0.0062?; 3060 + 1.72?; 
2718.3 - 1.52i 3885.6 + 5.0 x 10" 



(6) 



We need these matrices to calculate flavor changing processes. We find that all the fla- 
vor changing constraints arising from the SUSY exchange are consistent with the bounds 
obtained in [[TBI. The six down type squark masses in this example are given by: 



(459,458, 440, 439, 433 and 415) GeV and M- g = 501 GeV. The parameter space of the 
model is quite constrained, the example above and its overall rescaling (discussed later) are 
the only solution we have found. 

Let us now turn to CP violation in this model. Note that as mentioned before, the 
value of ex is used to determine the input phase of the theory. It is clear from the CKM 
mixing matrix of the above example that the rephasing invariant J-parameter is of order 
~ 7 x 10~ 7 . If the CKM phase is to explain ek, the value needed is J ~ 2 x 10~ 5 . Thus 
€k has a purely supersymmetric origin here. The dominant contribution is from gluino box 
graph involving the LL and RR terms in the squark mass matrix. These LL and RR terms 
also have their origin in the off-diagonal A terms through the RGE. All parameters in our 
model are then essentially fixed. In fact we have tried to vary them to see the effect. What 
we find is that fitting the quark mixings essentially implies that we must vary mo, min and 
A by a common factor k relative to the example just given. The value of ex is sensitive 
to k since the CKM contribution to the real part of Arrix is insensitive to it whereas the 
imaginary part of the matrix element which receives its dominant contribution from the 
supersymmetric box graphs scales like m^ 2 . Thus the only freedom allowed in our choice of 
squark and gaugino masses is whatever comes from the uncertainty in the hadronic matrix 
elements. Using the 12 elements of the LL and RR mass matrices of Eq. (6), we find 
that the experimental bound on ex is nearly saturated. We use Ref. [|17]] to find the QCD 



corrected bound on 



\Im(6 



! 2, LL ^12, RR) 



which is ~ 1.5 x 10 



for rriq 



460 GeV and 
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x ~ 1.2. (^ 2>iL = M? 2>LL /ml) 

Let us now turn to our prediction for e'/e. As usual the dominant contributions here 
are from the penguin diagrams. An important point to note is that since all parameters are 
now fixed, we might have a conflict with the measured value of e'/e. The reason for this 
apprehension is that Masiero and Silvastrini |T7j quote an upper limit on ImSu^LR — 1CT 5 , 
which is a factor of 10 smaller than the value required by us to fit ex- This is where 
exact parity symmetry comes to the rescue. The point is that the CP- violating AI = 1/2 
penguin Hamiltonian contributed by the exchange of squarks is proportional to the difference 
8n,LR — ^2i lr, which vanishes above the vr scale due to the constraints of parity symmetry 
mentioned above. However, the RGE running makes this difference nonzero by roughly 
~ j^ln Msusy ' wn ^ cn * s a ^ ^he l eve l- This qualitative conclusion is borne out by 
our detailed numerical calculation. The diagram for the operator d^a^t^s^G^ is formed 
by the squark line — hi, — sr in the loop formed by the squark and gluino lines. The 
magnitudes of the mixings are given <5i3,ll(= M^^/m-) and 8%2,lr{= /m~). Similarly 

we have the other diagram where the operator dRa^t^s^G^ is formed by the squark line 
dn — bR — sl- We have to subtract one diagram from the other to calculate e'/e and we 
predict e'/e ~ 3 x 10~ 3 for the above squark mass matrix and the lattice value for the 
hadronic matrix elements. 

Let us now discuss the parameter space where we can have the right amount of CP 
violation in e'/e. One simple way is if we change m , Mu 2 and the matrix A in our example 
by a common factor of k, then Vckm and the fermion masses will remain the same. However 
one needs to change the phase in order to fit €k but e'/e might go out of the experimental 
range. For example if we use k = 2, which corresponds to rrig = 900 GeV, we find that the 
phase is near 0.1 to fit ex, however our prediction of e'/e becomes a factor 2.5 smaller 
than what we had before. If k is decreased to 0.5 (corresponds to rriq = 230 GeV), then the 
prediction of e'/ e becomes a factor of 3 larger. The detailed predictions for the CP violating 
parameters e'/e along with the neutron edm d e n for various choices of k are given in Table I. 
We see from the Table that k somewhere between 0.7 to 2 is acceptable. This is the allowed 
spread in the squark mass parameters and the other SUSY breaking parameters i.e. squark 
and gluino masses somewhere between 300 GeV to a TeV. Note that the ratio of the gluino 
mass to the squark mass is essentially fixed, it is about 1.2 in all the fits. This prediction 
could serve as a crucial test of the model. We also find that tan/3 cannot be increased 
beyond about 6. Larger tanj3 would require larger value of the off diagonal elements of A 
(to fit V c b, V us etc). Through RGE, this would yield a SUSY contribution to Am K that is 
beyond the experimental limit. Note that the contribution to Amx will grow as (tan/3) 2 , 
so there is really very little room for tan (3 > 6. 

Coming to the electric dipole moment of the neutron, we first note that, at the v r scale, 
the flavor structure of the model is specified by the diagonal Yukawa matrices and Hermitean 
A matrices. The neutron edm would vanish in this limit, since it is given by the imaginary 
component of the (11) element of the A matrix. However once we extrapolate down to the 
weak scale, the situation changes and we get a nonvanishing, but quite consistent, edm for 
neutron. 
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Table 1. The predictions for — and neutron edm for different values of k. k 
corresponds to m = 80 GeV and mi/ 2 = 180GeV. 



k 




Neutron edm (ecm) 


0.5 


1.0 x 10~ 2 


5.0 x 10" 29 


1.0 


3.1 x 10~ 3 


5.2 x 10" 29 


1.5 


2.4 x 10~ 3 


6.0 x 10" 29 


2.0 


1.2 x 10~ 3 


7.0 x 10~ 29 



To calculate neutron edm we have considered 3 operators j4|: 0\ = — \qcf llv ^qF l 



2 = ~\qa^T a qGl u and 3 = |/ abc G» p G^G c A(T e^ A -, where G x yz is the gluon field strength 
and f a bc are the Gell-mann coefficients. The effective Lagrangian with the Wilson coefficients 
is given by: L = J2i=iCi(Q)Oi(Q). We evaluate the C^s at the weak scale and then we 
multiply by rji in order to evaluate them at 1.18 GeV ||. We use r\\ ~ 1.53 and r] 2 =r] 3 =3A. 
Finally we use naive dimensional analysis |nj to calculate the quark edms (d q = Ci(1.18) + 
^-C2(1.18) + ^^C3(1.18)) and then use the quark models to calculate the neutron edm. 
We estimate di 1 



10 —10 ecm. A rough intuitive way to see this number is to note 
that the dominant contribution to d e n comes from the An term which is complex and use 

the naive estimate from the formula |j£ {^^j Im[$ii,LR\- 

Before concluding, let us comment on other models with a supersymmetric origin of CP 
violation. One of the early models of this type is that of Ref. |JT9[] , where all CP violation is 
supposed to result from the phase of the gluino mass. Our model is different from that work 
since parity symmetry makes the gluino mass real above v r and small at the weak scale. A 
related class of models with approximate CP invariance [2(| generically leads to superweak 
CP violation which is now excluded. The second model is the recent phenomenological 
analysis of Ref. pip , where it is noted that a value of Su^lr — 10 -5 would produce the right 
value for the e'/e satisfying all other constraints. Here we have constructed essentially a 
complete theory that on RGE extrapolation leads to these parameters. As a result, we have 
more specific predictions e.g. the d e n as well as the fact that there will be no measurable 
CP violation in the B-sector. It should be noted that while the constraint of parity by 
itself does not require the Wino mass M 2 to be real, when embedded into a scheme such as 
SO(10) where M 2 gets related to M 3 , M 2 becomes real and there is no additional source of 
CP violation. 

In conclusion, we have found that the requirement that MSSM be embedded into a 
supersymmetric left -right framework above a scale of 10 12 GeV to explain the small neutrino 
masses observed in the Super-Kamiokande experiment, imposes very stringent constraints on 
the parameters of the MSSM. First it predicts the value of e'/e in agreement with experiment 
and the edm of neutron comfortably consistent with the present upper limits. We also find 
the super-partner masses to be in a very narrow range with the gluino mass not much 
different from the squark masses, which can provide a test of the model. The parameter 
tan j3 is not expected to exceed about 6. In the CP violation sector, a test of the model is the 
absence of significant CP violating effects in the B-sector. The standard KM contribution is 
negligible, so is the supersymmetric contribution (see Eq. (6)). These models can be readily 
unified into SO(10) models or other unification groups that contain the SUSY left-right 
model as a subgroup in them. 
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